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Abstract 


Agave species are high-yielding crassulacean acid metabolism (CAM) plants, some of which are grown commer- 
cially and recognized as potential bioenergy species for dry regions of the world. This study is the first field trial 
of Agave species for bioenergy in the United States, and was established to compare the production of Agave 
americana with the production of Agave tequilana and Agave fourcroydes, which are produced commercially in 
Mexico for tequila and fiber. The field trial included four experimental irrigation levels to test the response of 
biomass production to water inputs. After 3 years, annual production of healthy A. americana plants reached 
9.3 Mg dry mass ha! yr! (including pup mass) with 530 mm of annual water inputs, including both rainfall 
and irrigation. Yields in the most arid conditions tested (300 mm yr‘ water input) were 2.0- 
4.0 Mg dry mass ha! yr~'. Agave tequilana and Agave fourcroydes were severely damaged by cold in the first 
winter, and produced maximum yields of only 0.04 Mg ha! yr~' and 0.26 Mg ha“! yr |, respectively. The 
agave snout weevil (Scyphophorus acupunctatus) emerged as an important challenge for A. americana cropping, 
killing a greater number of plants in the higher irrigation treatments. Physiological differences in A. americana 
plants across irrigation treatments were most evident in the warmest season, with gas exchange beginning up to 
3 h earlier and water use efficiency declining in treatments with the greatest water input (780 mm yr ' water 
input). Yields were lower than previous projections for Agave species, but results from this study suggest that 
A. americana has potential as a bioenergy crop and would have substantially reduced irrigation requirements rel- 
ative to conventional crops in the southwestern USA. Challenges for pest management and harvesting must still 
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be addressed before an efficient production system that uses Agave can be realized. 
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Introduction 


Droughts limit agricultural production in arid and semi- 
arid regions, yet plants that exhibit crassulacean acid 
metabolism (CAM) can improve crop yields in these 
conditions (Cushman et al., 2015; Davis et al., 2015; Yang 
et al., 2015). Plants in the Agave (L.) genus are obligate 
CAM plants, some of which are produced commer- 
cially, and are being considered for bioenergy crops 
(Borland et al., 2009; Davis et al., 2011, 2014, 2015; Yang 
et al., 2015). In regions of the world where Agave is 
grown commercially for beverages or fiber, yields range 
from 8.5 to 22 Mg ha’ depending on species and loca- 
tion (Davis et al., 2014), but the latitudinal range for 
these crops is uncertain. Despite the potential advan- 
tages and high theoretical yields (Davis et al., 2011, 
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2014; Owen et al., 2015), there are many dry regions of 
the world including the southwestern USA where Agave 
has not yet been cultivated commercially or experimen- 
tally, making it difficult to assess the realistic potential 
of this genus for agricultural production. There are no 
field experiments to date that compare the yields of 
Agave species in arid to semiarid conditions or the 
response of Agave species to variable water inputs. 
Here, we provide the first field-scale yield estimates of 
Agave in the arid southwestern USA. 

Climate projections indicate that many semiarid and 
arid parts of the world will be subject to more frequent 
and extreme drought events in the future (IPCC 2014). 
Agricultural production in many of these areas is 
already vulnerable to drought, dependent on large vol- 
umes of irrigation, and susceptible to soil degradation 
(Scanlon et al., 2007). The yields of 8.5-22 Mg for Agave 
referenced above (Davis et al., 2014) reflect commercial 
production in Mexico and Tanzania where the climate 
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is semiarid due to prolonged periods of drought even 
though heavy rainfall does occur intermittently. As 
water use efficiency (WUE) is the primary advantage of 
CAM plants that has sparked interest in these plants for 
dryland agricultural production, there is a need to 
resolve how much precipitation would be needed to 
achieve commercially viable yields of Agave. Commer- 
cial production of bioenergy in the USA is primarily 
supported by Zea mays L. (corn) feedstock, which yields 
an average of 9.4 Mg ha! yr! on prime agricultural 
land, but the dominant commercial cropping systems in 
Arizona yield far less (eg, cotton yields 
1.5 Mg ha! yr’). 

The key physiological advantage of CAM, to distin- 
guish from C; and Cy, is the temporal separation of car- 
boxylase activities catalyzed by the enzyme PEPC 
(phosphoenolpyruvate carboxylase) and the enzyme 
RUBISCO (ribulose bisphosphate [RubP] carboxylase/ 
oxygenase) so that PEPC is primarily active during the 
night when stomata are open (opposite of C3 and Cy 
plants), and RUBISCO is primarily active during the 
day when light energy is available. With PEPC activa- 
tion occurring at night when temperatures are cooler, 
evaporative demand is lessened and the ratio of carbon 
dioxide to water that diffuses through the stomata is 
increased relative to what would be expected during 
daytime conditions. The carbon is stored as malic acid 
overnight in large vacuoles typical of succulent plants 
and then released by decarboxylation at dawn. This 
strategy for growth gives CAM plants a strong advan- 
tage over C, and C, plants in conditions of water stress 
(Osmond, 1978; Nobel, 1996; Winter & Smith, 1996; 
Davis et al., 2014). Many CAM plants also have thick 
waxy cuticles, sunken stomata, and roots that are 
responsive to water stress; all adaptations that hydrauli- 
cally isolate the plant in extremely dry conditions. Theo- 
retical projections indicate that high-yielding CAM 
plants such as Agave and Opuntia (Mill.) can contribute 
substantially to the global renewable energy portfolio 
(Davis et al., 2011, 2014; Cushman et al., 2015; Owen 
et al., 2015). 

There is growing interest in Agave species for bioen- 
ergy production in the arid and semiarid southwestern 
region of the USA. Commercial agricultural production 
has been demonstrated for species native to Mexico 
(e.g., Agave tequilana F.A.C. weber and Agave fourcroydes 
Lemair), but these have not been grown in the USA out- 
side of small horticultural and landscaping plantings. 
Agave americana L., also known as the American century 
plant, is native to the USA (native range includes Ari- 
zona, California, Florida, Louisiana, and Texas) and has 
been reported to reach sizes similar to its high-yielding 
Mexican relatives (2 m in height before flowering, >3 m 
in width) and achieve high photosynthetic rates under a 


range of environmental conditions (Neales, 1973; Gen- 
try, 1982; Nobel, 2003; Garcia-Moya et al., 2011), but 
field-scale yields of A. americana have not been assessed 
in previous literature. A. americana is produced for a 
beverage called pulque in Mexico (La Barre, 1938), but 
there is less information about the agricultural produc- 
tion of this species than what is available from the lar- 
ger-scale tequila and fiber industries. Tissue 
composition of A. americana has been measured, and the 
leaves of this species have a low percentage of lignin 
mass relative to other species of Agave as well as high- 
soluble carbohydrate concentrations relative to other 
advanced bioenergy crops (Mylsamy & Rajendran, 2010; 
Li et al., 2012; Corbin et al., 2015), both beneficial traits 
for biochemical conversion pathways to liquid fuel 
(Davis et al., 2011). 

The first US field trial of Agave species for bioenergy 
was established in 2012 in Maricopa, AZ, a desert agri- 
cultural area where annual rainfall averages about 
190 mm yr |. This study reports the results to date for 
this experiment. The agricultural plantings were 
designed as a fully factorial two-way experiment to 
determine differences in survivorship, productivity, and 
WUE among three species (A. americana, A. tequilana, 
and A. fourcroydes) and four irrigation treatments (rang- 
ing from arid to well-watered conditions), and to deter- 
mine whether responses to irrigation are species 
dependent. The study site is at the edge of the climate 
tolerance range for A. tequilana and A. fourcroydes 
(Lewis et al., 2015), but these are the large-scale com- 
mercial species against which a comparison of the US 
native species, A. americana, would be meaningful. Sig- 
nificant mortality and cold damage of A. tequilana and 
A. fourcroydes occurred after a rare cold event in the first 
year, so the focus of the study was narrowed by the 
third year to resolving differences in A. americana pro- 
duction, WUE, and survivorship across irrigation treat- 
ments. Early survivorship and production of 
A. tequilana, A. fourcroydes, and A. americana are 
reported, but differences in carbon dioxide assimilation, 
biomass yields, and WUE across four irrigation treat- 
ments were resolved only for A. americana in year 3. 


Materials and methods 


Study site and experimental design 


The study site is a 1.2-ha field located at the University of Ari- 
zona Maricopa Agricultural Center (MAC) in Maricopa, AZ, 
and is centered at 33°03/32" N latitude, 111°58’12” W longitude, 
and 361 m above sea level. The field was plowed, tilled, and 
leveled before planting. Soil in the field is classified as Casa 
Grande sandy clay loam, which is described as fine-loamy, 
mixed, superactive, hyperthermic Typic Natrargids (Post et al., 
1988). The field was subdivided into eight sections (each 
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approximately 32 m by 46 m) separated by berms approxi- 
mately 50 cm in height to contain water that was allocated to 
each treatment by way of flood irrigation (Fig. 1). Two repli- 
cates of four irrigation treatments were assigned randomly to 
the eight sections. Irrigation water was provided by a concrete- 
lined irrigation ditch located on the south end of the field. A 
2.5-m long piece of PVC pipe (305 mm, diameter) was hard- 
plumbed through a hole drilled in the irrigation ditch. Irriga- 
tion water to the field was controlled by a swivel valve located 
inside the ditch at the PVC pipe entry hole. The PVC pipe sys- 
tem was reduced to a 4.6-m long piece of 203-mm-diameter 
PVC pipe that contained an in-line propeller-driven flow meter. 
The PVC-valve system was then connected to a 203-mm-dia- 
meter plastic polypipe that ran along the center of the field in 
the north-south direction (Fig. 1). Gates installed along the 
polypipe at 1.0-m spacing allowed individual sections to be 
irrigated at desired times. Water was applied in approximately 
100-mm increments (the minimum needed to flood-irrigate the 
section uniformly), but exact amounts varied, as explained 
below. 

Irrigation treatments were initiated in the second year 
(2013), after planted sections had been uniformly irrigated in 
2012 (year one). The control treatment received only precipita- 
tion (200 mm annually on average from 2013 to 2014) and the 
minimal water needed to apply fertilizer. In addition to the 
native precipitation, the irrigation schedule and levels targeted 
for other plots were as follows: (i) twice annually (in March 
and July) to add ~200 mm yr7! (or ~2 MI ha“! yr7'), (ii) every 
other month from March to October to add ~400 mm yr! (or 
~4 MI ha™! yr7'), and (iii) monthly from March to October to 
add ~800 mm yr! (or ~8 Ml ha7! yr'). The actual water 
inputs achieved differed some from these target quantities due 
to variations in water flow rate, soil hydraulic and roughness 
conditions, integrity of berms, which occasionally required 


NORTH | | 


LOW IRRIGATION +100 mm 
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rebuilding after being broken by flood water, and practical 
decisions made by technicians during sporadic heavy rainfall 
events. Exact amounts of water applied during irrigation were 
measured (both flow rate and total volume) by the in-line flow 
meter. The actual 2-year (2013 and 2014) average amounts of 
irrigation applied for the four treatment levels were 
100 mm yr~!, 260 mm yr~!, 330 mm yr7!, and 580 mm yr7'. 
Including the mean annual rainfall (~200 mm; see detailed 
climate in Fig. S1), total water inputs in the treatments were 
300 mm yr, 460 mm yr, 530 mm yr, and 780 mm yr. 

The irrigation treatments described above were initiated in 
the spring of 2013 (year 2) and repeated through 2014 (year 3). 
The entire field was fertilized annually in June with N:P:K fer- 
tilizer in a ratio of 40:80:80 based on recommendations for com- 
mercial plantation management used in the tequila industry. 
Applied rates each year were 56 kg N ha ', 112 kg P ha’, 
and 112 kg K ha’. Weeding was accomplished manually in 
the first year, and required regular attention (twice per month). 
In the second and third years, weeds were controlled with her- 
bicide (Surflan) applied 3-4 time per year at a rate of 4.7— 
9.5 1 ha~', by mowing between rows where possible, and man- 
ual hoeing of weeds near the base of plants. 

Within each of the replicated irrigation sections, there were 
six plots of Agave for a total of 48 plots across the field (Fig. 1). 
Two plots each of A. tequilana, A. fourcroydes, and A. americana 
were randomly assigned within each replicated irrigation treat- 
ment, so that across the whole field there were four replicate 
plots of each species in each irrigation treatment. Each plot was 
an experimental unit and contained 49 plants (7 x 7) spaced 
2m apart both within and along rows. The experimental 
design in the field was purposeful to statistically test the differ- 
ence in survivorship, yield, and WUE among species and 
among irrigation treatments. We treated this as a fully factorial 
design instead of a split-plot design because each combination 


Fig. 1 Field layout of the experimental design with three species (AGAM: A. americana, AGFO: A. fourcroydes, and AGTE: A. tequi- 
lana) and four irrigation treatments. Each square plot contained 49 plants spaced 2 m apart. Four replicate plots of each species were 
randomly assigned to each irrigation treatment with two replicated and randomly located areas of the field assigned to each irrigation 


level (100, 260, 330, and 580 mm yr). Total water inputs including rainfall for each treatment were 300, 460, 530, and 780 mm yr. 


a 


The irrigation pipe (shown in gray) was laid through the center of the field to allow access to each section that would be flood-irri- 
gated at separate times. Berms (shown in brown) were constructed to separate each irrigated section. 
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of irrigation and species treatments was replicated throughout 
a single field and the sections delineated in the field for each 
irrigation were not unique. 

Agave tequilana plants were a single genotype of weber var. 
azul propagated clonally and purchased from growers in Oax- 
aca, Mexico. Agave fourcroydes plants were mixed unidentified 
genotypes purchased from growers in Yucatan, Mexico. Both 
A. tequilana and A. fourcroydes were imported bare rooted with 
permits from the U.S. Department of Agriculture and were uni- 
form in age of just under 1 year. Agave americana plants were 
from a variety of unidentified genotypes purchased from a 
grower in California, USA, and shipped to the site in one gal- 
lon pots (also approximately 1 year in age). 

All species were planted in the spring of 2012 between 
March 23 and May 5 and irrigated weekly for the first month 
to ensure establishment, then irrigated monthly until October 
2012 (~1000 mm water applied in total). Establishment was 
successful for all three species, and there was no mortality 
observed prior to January of 2013. There was an unusually cold 
period with nine consecutive days of temperatures below freez- 
ing in January of 2013, and minimums between —6 °C and 
—8 °C for five consecutive nights. The A. tequilana and A. four- 
croydes plots suffered severe tissue damage. The monthly mean 
minimum temperature was 0.6 °C, 67% lower than the mean 
minimum for January in the other growth years. Although 
many plants remained, they were severely damaged and the 
biomass of both A. fourcroydes and A. tequilana were magni- 
tudes lower than the biomass of A. americana as a result. The 
study is therefore focused on the response of A. americana to 
irrigation, although survivorship and yields were quantified 
and compared among species after year two. There were an 
insufficient number of plots of A. tequilana remaining across 
irrigation treatments to statistically assess the response to irri- 
gation for this species. 


Measurements 


A field-calibrated neutron probe (Model 503, Campbell Pacific 
Nuclear, CPN, Martinez, CA) was used to measure the volu- 
metric soil water contents (0) of plots from 0.1 m to 1.9 m in 
0.2 m incremental depths. Between mid-April and mid-May, 
2013, 33 neutron access tubes were installed vertically to a soil 
depth of 2.0 m using a tractor-mounted Giddings soil sampler 
(Model 25-TS; Giddings Machine Company, Windsor, CO, 
USA). Of the 33 tubes, 16 were located within each of the 16 
A. americana plots (i.e., four per irrigation level). Tubes were 
installed at a distance of about 1.0 m away from plants. During 
installation of the neutron access tubes, soil samples in 0.3-m 
increments were collected at each location to a depth of 1.8 m 
(with 6 depth increments). Each separate soil sample was ana- 
lyzed for soil particle size fraction using the Bouyoucos 
hydrometer method (Gee & Bauder, 1986). Although sandy 
clay loam and sandy loam textures were predominate, soil tex- 
ture analyses revealed approximately 10% of the total samples 
were classified either as clay loam, loam, or sandy clay. The 
volumetric soil water content measurements at the 33 access 
tube locations began on May 21, 2013. In general, soil water 
contents were measured twice monthly, just before and a few 


days after a scheduled irrigation. Soil water contents were mea- 
sured for all access tubes and depths with each irrigation event, 
even when only a subset of the plots were irrigated. 

Evapotranspiration (ET..) for the irrigation treatments of 
A. americana was determined over successive soil water content 
measurement dates as the residual of the soil water balance of 
an estimated root zone of 2.0 m. Following Jensen et al. (1990), 
the soil water balance for ET., occurring between the first mea- 
surement date, denoted as day i = 1, and the second date ~ 
15 days later on day n was calculated as: 


10 n—1 


ET, = X (Si — Sin) + DR at IW;) (1) 


i=l j=l 


where ET.. is the total evapotranspiration that occurred from 
day 1 to day n, S;ı and S; „ are soil water storage measured at 
each of the 10 neutron probe soil depth layers on day i = 1 and 
on day n, respectively. R; and IW; are, respectively, the mea- 
sured rainfall and applied irrigation depths received on day j. 
All units are in mm. Rainfall data were provided by a University 
of Arizona, Meteorological Network (AzMet; ag.arizona.edu/ 
azmet) weather station located at MAC, about 1.2 km from the 
field site. For each of the 10 soil depth increments, soil water 
storage was calculated as the measured soil water contents (m?/ 
mî) at the layer times the soil depth layer in m (i.e., 0 x 0.2 m). 
For each soil water measurement location, the volumetric soil 
water contents at field capacity (FC) and permanent wilting 
point (PWP) were determined for each layer from the soil tex- 
ture data using soil water characteristics estimation procedures 
developed by Saxton & Rawls (2006). Equation 1 does not 
include a deep percolation component, as there was no evidence 
of soil water content increasing below a soil depth of 1.5 m. 
Fig. 2 shows the periodically measured 0 for each layer with 
time for the wettest irrigation treatment (i.e., 580 mm). The fig- 
ure includes the FC and PWP for each layer determined from 
soil texture. For 2013 (Fig. 2a), soil water contents measured a 
few days after irrigations increased to a depth of 1.5 m. How- 
ever, the soil water contents at lower layers did not increase 
above FC, indicating negligible drainage from the root zone. 
Conversely, it can be seen, particularly in 2014 (Fig. 2b), that soil 
water was depleted somewhat to 1.9 m, suggesting that a root 
zone of 2.0 m for agave water use was reasonable. 

Survival and reproduction rates of plants were recorded and 
compared across species and irrigation treatments in the second 
year of growth. Survival rates were recorded in July of 2013 
(after all cold damage could be assessed) and biomass produc- 
tivity was recorded in January of 2014 (after two full growing 
seasons). The mean percentage of surviving plants in four repli- 
cate plots was used to estimate survivorship of each species in 
each irrigation treatment. No significant cold damage was 
observed after the spring of 2013. Survivorship and biomass 
productivity was measured again in February of 2015 (after 
three full growing seasons). The mean number and condition of 
reproductive offsets (pups) was also measured in each plot. 

The response of dusk-to-dawn carbon dioxide uptake to the 
different irrigation treatments was determined by measuring 
gas exchange in A. americana using Li-Cor LI-6400XT Portable 
Photosynthesis Systems (Li-Cor; Lincoln, NE, USA). These 
infrared gas analyzers were used to record carbon uptake rates 
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Fig. 2 Soil water content measurements from 0.1 to 1.9 m, in 
0.2-m increments (dotted lines), field capacity (FC, solid line), 
and permanent wilting point (PWP, dashed line) derived from 
soil texture analyses for 0.3-m soil layers from 0 to 1.8 m for 
the 800-mm treatment plot in 2013 (a) and 2014 (b). Only a sub- 
set of all measurements are shown as examples. Note soil water 
content measurements made after an irrigation are denoted by 
dark symbols. The FC and PWP data are plotted at the center 
of the layer depth (e.g., at 0.15 m for the 0- to 30-m layer). 


for plants in each irrigation treatment over a time course start- 
ing before sunset and extending until after dawn at three times 
throughout the year (March, June, and November). Total net 
carbon assimilation was estimated by integrating the measure- 
ments of gas exchange over time (typically a 15-h period). Sam- 
pling began in the late afternoon in attempt to include gas 
exchange during Phase IV, but gas exchange was not observed 
prior to 6 pm. In November of 2014, a full 24-h time period 
was measured to verify that the 15-h measurements reflected 
the full duration of stomatal opening. Dynamics of gas 
exchange were also compared to determine whether the timing 
of peak carbon uptake rates and the phases of CAM differed 
across the four different irrigation treatments. 

During each seasonal sampling, gas exchange over the 15-h 
period from 6:00 pm to 9:00 am was measured at approxi- 
mately 30-min intervals in each irrigation treatment. Each 
incremental estimate during the diurnal time period was a 


mean of the gas exchange rate measured in a mature leaf 
located midway between the outside whorl and center spike of 
three randomly selected plants in each irrigation treatment. 
This approach allowed us to estimate gas exchange levels using 
a mean of several plants within an irrigation treatment instead 
of that of a single plant. A square 2 cm x 3 cm cuvette was 
used with thick gaskets to accommodate the thick succulent 
leaves, and measurements were made at a point on the leaf 
approximately two thirds the distance from the base. There 
were occasional gaps in data (a missing 30-min increment) due 
to the need to switch batteries in the LiCor 6400 or other main- 
tenance requirements. When calculating the total dusk-to-dawn 
CO, assimilation, measurements from the next closest time 
interval were used to estimate fluxes during each gap. 

Aboveground biomass of a random plant from each of the 
sixteen A. americana plots was measured after destructively 
harvesting at the end of second and third years of growth 
(mid-January 2014 and early February 2015). In February of 
2015, aboveground biomass of all plants that were killed by the 
agave snout weevil was also measured (the number of plants 
killed per plot ranged from 3 to 38). Because most plants har- 
vested in February of 2015 were infested with the snout weevil, 
an additional harvest occurred in the spring of 2015 to assess 
the productivity of plants that remained healthy. Subsamples 
from each plant, consisting of three randomly selected leaves, 
were dried in an oven for 2 weeks at 65 °C and then weighed 
to estimate plant tissue water content. The water mass was sub- 
tracted from the total wet biomass to estimate dry biomass pro- 
duction in each plot. 

Net WUE was estimated in two ways: WUE), as total dry 
biomass per unit of water input (including both precipitation 
and irrigation); and WUE, as total dry biomass per unit ET.. as 
calculated from the soil water balance described above. WUE; 
provides a practical estimate for growers of the amount of bio- 
mass one can expect given a certain water input. WUE, can be 
used to resolve any additional soil water that might be used to 
supplement the plant water requirements. 

The statistical program JMP Pro 10 (SAS Institute Inc., Cary, 
NC, USA) was used to assess the differences in survivorship 
and biomass productivity in year 2 among species and irriga- 
tion treatments, and whether the effect of irrigation depended 
on species, using a two-way analysis of variance (anova). The 
same statistical program was used to test differences in (i) sur- 
vivorship, (ii) reproductive offsets, (iii) biomass productivity, 
and (iv) WUE of A. americana among irrigation treatments in 
year 3 using a one-way Anova. The Shapiro—Wilk W-test was 
used to determine whether the distribution in response vari- 
ables met the assumption of normality. Tukey’s HSD (honest 
significant differences) was used to resolve significant differ- 
ences among treatment levels. Linear regression was used to 
compare the two estimates of WUE, and nonlinear regression 
was used to assess the response of WUE to increasing water 
inputs and ET.. 


Results 


Initial establishment of A. americana, A. tequilana, and 
A. fourcroydes was successful, with a mean survival at 
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the plot level of 100%, 97%, and 92%, respectively, by 
August of 2012. The few plants that were lost during 
this period were replaced. In January of 2013, the mini- 
mum temperature on five consecutive nights ranged 
from —6 °C to —8 °C, and cold damage was immedi- 
ately evident in the A. tequilana and A. fourcroydes even 
though mortality was not entirely evident until March 
of 2013. A survey of the field in July of 2013 revealed 
that 74% of the A. tequilana plants were killed by the 
prolonged cold (Fig. 3), and 88% of the A. fourcroydes 
were severely damaged. There was no mortality or cold 
damage observed in the plots of A. americana after the 
cold event in January of 2013, and survival of the 
A. americana remained near 100% until 2014. There was 
greater mortality in A. tequilana than in A. fourcroydes or 
A. americana, and irrigation treatments did not affect 
survival in the second year for any species. Statistical 
significance is not reported for these differences because 
the distribution of response variables was bimodal (vio- 
lating the assumption of normality in the statistical 
model) and could not be corrected through transforma- 
tion. Ranking data for analysis using nonparametric 
tests was not possible because survival was 100% for all 
A. americana plots (insufficient range for ranking). Nev- 
ertheless, the difference in survivorship among species 
was Clear (Fig. 3). 

Although many of the A. fourcroydes survived, the 
cold damage was severe and the aboveground tissue 
completely died back in many of the plots. After recov- 
ery, the aboveground biomass of A. fourcroydes in year 
two was less than 0.3 Mg ha! yr) even in the highest 
irrigation treatment, and was less than 9% of the bio- 
mass production measured for A. americana in all irriga- 
tion treatments (Fig. 4). While there was a significant 
main effect of irrigation on biomass (P < 0.01), this 
response differed by species (P < 0.01), with A. ameri- 
cana biomass increasing more with irrigation levels than 
biomass of A. fourcroydes. A. americana biomass was also 
significantly greater than A. fourcroydes biomass in all 
irrigation treatments (P < 0.01). Before statistically 
resolving the differences between A. americana and 
A. fourcroydes, biomass data were log-transformed to 
achieve a normal distribution. There were so few A. te- 
quilana surviving that replicate plots of living plants no 
longer remained in all treatments (some plots had 0 sur- 
viving plants), but the biomass production of the 
remaining plots was <1% of A. americana production 
(Fig. 4). A. tequilana and A. fourcroydes did not recover 
nor make any significant gain in biomass production by 
the end of the third year. 

The first evidence of mortality in the A. americana 
plots was recorded in January of 2014, and this was due 
to an infestation by Scyphophorus acupunctatus (agave 
snout weevil). The plots were first treated with an 


insecticide (Merit, Bayer CropScience, Monheim, Ger- 
many) in February of 2014, and two repeat applications 
occurred in the summer and fall (each at a rate of 67— 
87 kg ha” !). The weevil was, however, never fully erad- 
icated from the field. In February of 2015, 77% of the 
originally established plants (mother plants) in the high- 
est irrigation treatment had been killed by the snout 
weevil, significantly more than the other irrigation treat- 
ments (P < 0.01). Survivorship was greatest (80%) in the 
plots with the lowest annual water inputs (Fig. 5). 
Asexual reproduction rates of A. americana were 
greater in the plots with 780 mm annual water inputs in 
2014, but the mean number of reproductive offsets 
(pups) produced for each original plant (mother plant) 
was similar across treatments in 2015 (Fig. 6). An aver- 
age of 3 pups per mother plant were produced across 
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Survivorship in year 2 (%) 


A. americana A. fourcroydes A. tequilana 


Fig. 3 Percent survivorship of three Agave species by July of 
2013. Bars represent means of replicate plots across all irriga- 
tion levels (n = 16), and error bars represent standard error. 
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Fig. 4 Annualized productivity of three Agave species across 
the four irrigation treatments by January of 2014. Bars repre- 
sent treatment means with plots as the replicate (n = 4), and 
error bars represent standard error. Different letters indicate 
statistically significant differences (P < 0.05). 
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all plots of A. americana in 2014, and a mean of 4 pups 
per mother plant were produced by 2015. There was not 
a significant difference among irrigation treatments in 
the number of surviving pups remaining in the plots 
(P = 0.14), or in the total number of pups produced per 
mother plant (P = 0.10), although the trend in 2015 was 
toward greater pup production with annual water 
inputs of 460 mm yr! and 530 mm yr. Reproductive 
offsets were not observed in the other species. 

Biomass of A. americana plants increased with water 
inputs of up to 530 mm in all years (Fig. 7). In 2014, 
plots with the greatest water inputs (530 mm yr ' and 
780 mm yr‘) resulted in the highest annual biomass 


100 


A. americana survivorhip (%) 


300 460 530 780 
Irrigation + rainfall (mm yr-*) 


Fig. 5 Percent survivorship of A. americana after the agave 
snout weevil spread through the field in 2014. Survivorship 
was assessed in February of 2015. Bars represent treatment 
means with plots as the replicate (n = 4) and error bars repre- 
sent standard error. * indicates statistically significant differ- 
ence (P < 0.05). 


2013 m2014 2015 


Pups per mother plant (#) 


300 460 530 780 
Irrigation + rainfall (mm yr-*) 


Fig.6 Number of asexually propagated offsets (pups) 
observed for each A. americana plant that was initially estab- 
lished in the field (mother plants). Bars represent mean of repli- 
cate plots (n = 4), and error bars represent standard error. 


production (2.8-3.0 Mg ha ! yr’, significantly greater 
than the 300-mm or 460-mm treatments (P < 0.01). The 
difference among treatments in biomass of dead plants 
harvested in 2015 was not significant (P = 0.77). There 
was, however, a trend for greater biomass productivity 
of healthy mother plants in the plots that received 
530 mm yr’ of water when compared to those in the 
other treatments (P = 0.06; Fig. 7). Total aboveground 
biomass, including biomass of both pups and mother 
plants, was significantly different across irrigation treat- 
ments in 2015 (P = 0.02), with the biomass production 
reaching 9.3 Mg ha ' yr ' in plots with 530 mm of 
annual water input (Fig. 7). After 3 years, annual pro- 
duction of healthy plants was more than double the 
annual production after 2 years, although mean produc- 
tion rates of plants killed by the snout weevil were 
much lower (Fig. 7). 

WUE;, calculated as total aboveground biomass pro- 
duction per unit of total water input, was similar in irri- 
gation treatments that ranged from total inputs of 


10 E Mother plants 2014 (a) 
m Mother plants harvested dead 2015 
® Mother plants harvested live 2015 
® Total aboveground production of health 


Annual dry biomass yield (Mg ha~? yr-?) 


Dry biomass yield (Mg ha~?) 


300 460 530 780 
Irrigation + rainfall (mm) 


Fig. 7 Agave americana oven-dried biomass production annual- 
ized (a) and total (b) for 2-year-old mother plants (blue, 2014), 
dead 3-year-old mother plants (orange, 2015), living 3-year-old 
mother plants (gray, 2015), and total aboveground biomass of 
healthy mother and pups (yellow, 2015). Bars represent a mean 
of individuals from replicate plots (n = 4), but the number of 
dead plants per plot varied (3 < n < 31) and error bars repre- 
sent standard error. 
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300 mm yr’ to 530 mm yr! (P > 0.05), but WUE: 
declined in the highest irrigation treatment and was sig- 
nificantly lower with inputs of 780 mm yr ' (P = 0.01; 
Fig. 8a). Trends in the response of WUE», calculated as 
total aboveground biomass production per unit of ET.. 
(Fig. 8b), to water inputs were very similar to the 
response of WUE; because there was little ET., derived 
from soil water storage (Table 1). Results from the two 
methods estimating WUE were closely correlated 
(Fig. 8d; R? = 0.99). There was a nonlinear correlation 
between water inputs and WUE, (R? = 0.88) and 
between ET.. and WUF, (R? = 0.82), with the greatest 
WUE resulting from conditions with total water inputs 
of 460 mm yr' to 530 mm yr’ (Fig. 8). 

The greatest difference in carbon assimilation among 
treatments, as measured by dusk-to-dawn gas exchange, 
occurred in June (Fig. 9), which was the time period 
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with the lowest native rainfall (relative to March and 
November). During the dry summer nights, plants 
growing in plots with higher water inputs (530 mm and 
780 mm) opened stomata up to 3 h earlier than plants 
growing in the control plots (Fig. 10). Carbon assimila- 
tion increased with irrigation treatment levels in June, 
but differences across treatments were not evident in 
March and November (Fig. 9). The measurement cam- 
paign in November confirmed that there was no mid- 
day carbon assimilation. 


Discussion 


Agave americana was successfully established as an agri- 
cultural planting in the arid conditions of the south- 
western USA. Although biomass production of 
A. americana responded to irrigation, there was no bene- 
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Fig. 8 Relationships between increasing annual water inputs and water use efficiency calculated as the biomass per unit of water 


inputs, WUE; (a); between increasing annual crop evapotranspiration calculated from soil water balances, ET., and water use effi- 


ciency calculated as biomass per unit of ET, WUE; (b); and between increasing annual water inputs and oven-dried biomass (c) of 


A. americana. Each point in panels (a), (b), and (c) represents a mean of individuals from replicate plots (n = 4), open circles indicate 
biomass or WUE of mother plants only, and closed circles indicate biomass or WUE of all aboveground plant mass including mother 
plants and pups. Data here are representing only healthy living plants. Also shown is the close relationship between annual water 


input (rainfall + precipitation) and annual ET, (d). Each point in panel (d) represents a mean of replicate irrigation sections in the 


experimental field (N = 2). All error bars indicate standard error. 
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Table 1 Treatment means* for Agave americana soil water bal- 
ance components} expressed as yearly averages for 2013-2014, 
where IW is exact irrigation applied without rounding, R is 
measured rainfall, AS} is the change in soil water storage 
within a soil profile of 2.0 m, and ET. is annual crop evapo- 
transpiration 


Treatment IW (mm) R (mm) AS (mm) ET. (mm) 
580 577 202 14 793 
330 333 202 8 543 
260 261 202 36 499 
100 96 202 34 332 


*Treatments are denoted as amount of annual irrigation 
applied. 

+Note deep percolation water penetrating below 2.0 m was 
assumed negligible for all treatments. 

{AS is the change averaged from Jan. 1 to Dec. 31 in 2013 and 
from Jan. 1 to Dec. 31 in 2014. 


fit of irrigating beyond 330 mm yr! and the optimum 
water input (including rainfall) for an A. americana crop 
was ~530 mm yr. Yields under these conditions ran- 
ged from 2.5 to 9.3 Mg ha`' yr! (dry mass) over 
3 years depending on pest infestations and whether 
harvests included only mother plants or all plant bio- 
mass. The lower production estimate is comparable to 
that of soybean (~2.9 Mg ha~ yr), a crop grown for 
biodiesel, and the upper production estimate is compa- 
rable to corn grain (~9.4 Mg ha” yr~'), a feedstock for 
bioethanol production and the most abundant biofuel 
crop globally, in prime agricultural regions of the USA 
based on the national average of these crops over the 
last 5 years (NASS, 2015). 

In the most arid conditions tested here (300 mm of 
annual water inputs), biomass production of A. ameri- 
cana ranged from 2.0 to 4.0 Mg ha’! yr‘, indicating 
that A. americana crops may have commercially viable 
yields in dry regions with minimal irrigation required 
for establishment. Approximately 30% of global land is 
arid or semiarid, most of which is considered poorly 
suited or unsuitable for agriculture (van Velthuizen 
et al., 2007), and Agave crops offer opportunity for pro- 
duction in these regions. Irrigation for agriculture con- 
sumes more water than any other practice and 
continued high consumption is unsustainable, especially 
in dry regions (Giovannucci et al., 2012). 

In the southwestern USA, agriculture that relies heav- 
ily on irrigation is economically unstable. Cotton, for 
example, one of the most important crops in Arizona, 
requires an average of 1,046 mm yr' of irrigation 
water (in addition to rainfall) and yields only 
1.46 Mg ha ! annually (USDA 2012). This high water 
demand is increasingly problematic; recent estimates 
suggest that 20-60 Mha of irrigated cropland globally 
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Fig. 9 Net CO, uptake of A. americana measured at approxi- 
mately 30-min intervals between the hours of 6 pm and 9 am 
in March (light gray), June (green), and November (dark gray) 
of 2014. Each bar represents a mean of individuals (n = 3 for 
most time points) in each irrigation treatment, and error bars 
represent the standard error integrated over the full 15 h sam- 
pling period. Note that 29 mm of rain occurred prior to the 
March measurements, soil conditions were dry during the June 
measurements, and the agave snout weevil was disproportion- 
ately present in the higher irrigation treatments in November. 


will be forced out of production by the end of the cen- 
tury due to a decline in the availability of freshwater 
resources (Elliott et al., 2014). A. americana is an alterna- 
tive crop that can be cultivated in dry climates with as 
little as a tenth of the water inputs required from con- 
ventional crops. 

Although moderately productive with 300 mm - 
460 mm of annual water input, the soil water balance 
measurements (Table 1) suggest that a greater amount 
of stored soil water was withdrawn by A. americana in 
these conditions (Fig. 8d). The additional water demand 
is small (34-36 mm), but could have an impact on soil 
water recharge rates in dry regions if commercially 
scaled plantations were established. Effects of agricul- 
ture on groundwater recharge are a major challenge for 
the southwestern USA (e.g., Scanlon et al., 2005). The 
optimum annual water input therefore seems to be 
between 460 mm and 530 mm of total water inputs. At 
inputs >530 mm, ET.. is nearly equal to water inputs 
(Fig. 8d). 

The agave snout weevil is a major threat to crops of 
A. americana, and the impact of the snout weevil on 
A. americana survival increased with irrigation levels. 
The susceptibility of the plants in the high irrigation 
treatments could be due to the larger initial size of these 
plants. Horticultural texts indicate that the weevils pref- 
erentially burrow into larger plants to reproduce (Irish 
& Irish, 2000). Alternatively, and because the plant mor- 
tality due to the snout weevil was observed in plants 
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that were only 2 years old, the additional moisture in 
the high irrigation treatments may provide an optimum 
environment for the agave snout weevil. 

The agave snout weevil is a well-known threat to 
A. americana plants, but effective treatments have not 
been fully developed. Recent studies indicate that 
extracts from castor oil plants (Ricinus communis L) may 
be effective repellents for adult weevils (Pacheco- 
Sanchez et al., 2012a,b). Adult weevils burrow into the 
stem of the Agave plant and reproduce there, so the lar- 
vae are nourished by the readily available carbohy- 
drates at the leaf bases. The result is that the base of the 
stem is destroyed before the rest of plant shows symp- 
toms (Fig. 11). The agave snout weevil is widespread in 
the United States, and should be expected in an A. amer- 
icana crop unless a resistant genotype is identified. 
There is some anecdotal evidence of resistant genotypes 
(Irish & Irish, 2000), and it is possible that resistant indi- 
viduals will be identified as this field trial continues. 

If the timing of the snout weevil infestation can be 
delayed until at least the third or fourth growing sea- 
son, which seems most likely in conditions with low 
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water inputs, biomass yields can still be competitive 
with other bioenergy crops. It is noteworthy that the 
plants are easier to harvest manually after the snout 
weevil has weakened the stems, but minimizing the 
spread of this pest would be most desirable. Commer- 
cial production in the USA will likely require mecha- 
nized harvesting instead of the traditional manual 
harvesting that is common in Mexico. A recent eco- 
nomic analysis indicated that manual harvests were too 
costly for economically viable production of Agave for 
biofuel in the USA (Nuñez et al., 2011). Existing mecha- 
nized harvesters could be modified to cut the plants at 
the base as in tequila plantations or to trim leaves annu- 
ally as in fiber plantation. 

This study is the first field experiment to compare 
Agave production under a range of controlled water 
inputs. The production of Agave species measured in 
this study was lower than previously reported Agave 
biomass in semiarid regions (e.g., Davis et al., 2011), but 
A. tequilana production was similar to modeled projec- 
tions of A. tequilana biomass potential in the United 
States (Owen et al., 2015). There are no previous reports 
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Fig. 10 Dusk-to-dawn CO, uptake of A. americana measured in June of 2014 with four levels of water inputs (290 mm yr", yellow; 
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410 mm yr~', blue; 520 mm yr~', orange; 680 mm yr7!, gray). 


Fig. 11 Stem of A. americana affected by the agave snout weevil (left) and the stem of a healthy A. americana (right). 
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of A. americana production, but yields were lower than 
theoretical production rates for CAM plants (Nobel, 
1991; Borland et al., 2009; Davis et al., 2014), and far 
below the theoretical predictions of 43 Mg ha! yr for 
high-yielding CAM plants (Borland et al., 2009). Davis 
et al. (2014) projected that, in arid conditions with only 
200 mm annual precipitation, optimum yields of CAM 
plants could be 8.9 Mg ha! annually. Although the 
lowest irrigation treatment in this study, with 300 mm 
total water input, yielded 2.0-4.0 Mg ha! yr, it is 
possible that a more mature crop (24 years), optimized 
genetic lines of A. americana, or other cold-tolerant 
Agave spp. will yield greater biomass. 

Changes in planting density can have substantial 
impact on biomass production per unit area. The plant- 
ing density of each plot in this study was ~2500 plants/ 
ha. This is consistent with plant spacing used in sisal 
production and tequila plantations, and allows enough 
space for field managers to move between rows. If 
planting density could be increased to 4000 plants/ha, 
and managed with specialized equipment or remote 
technology, then higher yields could be expected. The 
crop canopy was not closed (there was exposed ground) 
in all plots measured in this study. 

The results of this experiment to date indicate that 
A. tequilana and A. fourcroydes, two important commer- 
cial species in Mexico, were not productive 3 years after 
planting in Maricopa, AZ. This region has been pre- 
dicted to have marginally appropriate growing condi- 
tions for A. tequilana, while more optimal conditions are 
predicted in Texas, California, and Florida (Lewis et al., 
2015; Owen et al., 2015). Still, the sensitivity of these spe- 
cies to cold make them vulnerable to sporadic extreme 
climate events, and additional risk would be introduced 
when planting these species in the southwestern USA. 
Variable cold tolerances among Agave spp. have been 
observed and described in previous literature (e.g., Nobel 
& Smith, 1983). The response of the tonoplast to tempera- 
ture increases during the day has been implicated in con- 
trolling phase changes that are essential to CAM (Kluge 
& Schomburg, 1996), but research to resolve mechanisms 
of cold tolerance in CAM plants is needed. Carbohydrate 
concentrations, particularly high molecular weight fruc- 
tans, are correlated with freeze tolerance in grasses 
(Dionne et al., 2010), and may explain the difference in 
cold tolerance between Agave species in this study 
because A. tequilana leaves have lower concentrations of 
water-soluble carbohydrates than A. americana leaves (Li 
et al., 2012; Corbin et al., 2015). 

Past life-cycle assessment indicated that there was 
potential for Agave to have equal or greater environmen- 
tal benefits for displacing fossil fuels when compared to 
sugarcane (Saccharum officinarum L.), corn (Zea mays L.), 
and switchgrass (Panicum virgatum L.) (Yan et al., 2011). 


Yan et al. (2011) assumed production rates of approxi- 
mately 30 Mg ha’ after 6 years (5 Mg ha! yr`') in 
A. tequilana stems, the sugar content of which was 
assumed to be 80% of the dry mass. Corbin et al. (2015) 
recently estimated the sugar content of A. americana 
plant tissue to be 60% of dry mass, suggesting that 
bioenergy production from A. americana in the south- 
western USA could have similar life-cycle benefits if the 
optimum yields (6.1-9.3 Mg ha! yr~') observed in this 
study were achieved. The feasibility of a full commer- 
cial production system that uses A. americana as feed- 
stock depends on other economic factors specific to the 
location where it would be cultivated. 


Conclusions 


Agave americana is a potential alternative bioenergy crop 
for the southwestern USA, even in arid conditions. A. te- 
quilana and A. fourcroydes were vulnerable to extreme 
cold events that can occur in this region, and were far 
less productive than A. americana in this side-by-side 
trial. The agave snout weevil poses a clear risk to the pro- 
duction of A. americana, but it may be manageable in 
conditions where annual water inputs are <530 mm. 
Genotypic selection for resistance to the agave snout 
weevil might also reduce this problem in the future. 
Healthy 3-year-old A. americana plots yielded 4.0— 
9.3 Mg ha! yr! total biomass with 300-530 mm yr! 
of water input (including 100-330 mm yr! irrigation, 
respectively), and plot yields without pup mass included 
were 3.1-6.1 Mg ha! yr. The production rates of 
A. americana were greater than conventional crops in Ari- 
zona (e.g., cotton) with less water input than is typical in 
this dry agricultural region. 


Acknowledgements 


We thank Glenda Simer, Richard Simer, Greg Main, and the 
Maricopa Agricultural Center for managing the experimental 
field. We also thank Michael Kampwerth and Michael Masters 
for helping with the initial establishment of the experimental 
field and appreciate assistance from Jose Ignacio Del Real 
Laborde in obtaining seedlings for the experiment. This project 
was funded in part by the Energy Biosciences Institute, and we 
appreciate valuable support and feedback from Steve Long and 
Chris Somerville. 


References 


Borland AM, Griffiths H, Hartwell J, Smith JAC (2009) Exploiting the potential of 
plants with crassulacean acid metabolism for bioenergy production on marginal 
lands. Journal of Experimental Botany, 60, 2879-2896. 

Corbin K, Byrt CS, Bauer S et al. (2015) Prospecting for energy-rich renewable raw 
materials: Agave leaf case study. PlosOne, 10, e0135382. doi: 10.1371/jour- 
nal.pone.0135382 

Cushman J, Davis SC, Yang X, Borland A (2015) Development and use of bioenergy 
feedstocks for semi-arid and arid lands. Journal of Experimental Botany, 66, 4177— 
4193. 


© 2015 The Authors. Global Change Biology Bioenergy Published by John Wiley & Sons Ltd., doi: 10.1111/gcbb.12324 


12 S.C. DAVIS et al. 


Davis SC, Frank G, Long S (2011) The global potential for Agave as a biofuel feed- 
stock. GCB Bioenergy, 3, 68-78. 

Davis SC, LeBauer D, Long S (2014) Light to liquid fuel: theoretical and realized 
energy conversion efficiency of plants using Crassulacean Acid Metabolism 
(CAM) in arid conditions. Journal of Experimental Botany, 65, 3471-3478. 

Davis SC, Ming R, LeBauer D, Long S (2015) Toward systems-level analysis of agri- 
cultural production from crassulacean acid metabolism (CAM): scaling from cell 
to commercial production. New Phytologist, 208, 66-72. 

Dionne J, Rochefort S, Huff DR, Desjardins Y, Bertrand A, Castonguay Y (2010) Vari- 
ability for freezing tolerance among 42 ecotypes of green-type annual bluegrass. 
Crop Science, 50, 321-336. 

Elliott J, Derying D, Miller C et al. (2014) Constraints and potentials of future irriga- 


tion water availability on agricultural production under climate change. PNAS, 
111, 3239-3244. 

Garcia-Moya E, Romero-Manzanares A, Nobel PS (2011) Highlights for agave pro- 
ductivity. GCB Bioenergy, 3, 4-14. 

Gee GW, Bauder JW (1986) Particle size analysis. In: Methods of Soil Analysis, Part 1, 
Physical and Mineralogical Methods, Agronomy Monograph No 9, 2nd edn (ed. 
Klute A), pp. 383-411. American Society of Agronomy, Madison. 

Gentry H (1982) Agaves of Continental North America. The University of Arizona 
Press, Tucson, AZ. 

Giovannucci D, Scherr SJ, Nierenberg D, Hebebrand C, Shapiro J, Milder J, Wheeler 
K. (2012) Food and Agriculture: The Future of Sustainability (March 1, 2012). The Sus- 
tainable Development in the 21st Century (SD21) Report for Rio+20. United Nations, 
New York. http://dx.doi.org/10.2139/ssrn.2054838 

IPCC (2014) Climate Change 2014: Impacts, Adaptation, and Vulnerability. Part A: Global 
and Sectoral Aspects. Contribution of Working Group II to the Fifth Assessment Report 
of the Intergovernmental Panel on Climate Change (eds Field CB, Barros VR, Dokken 
DJ, Mach KJ, Mastrandrea MD, Bilir TE, Chatterjee M, Ebi KL, Estrada YO, Gen- 
ova RC, Girma B, Kissel ES, Levy AN, MacCracken S, Mastrandrea PR, White 
LL). Cambridge University Press, Cambridge, United Kingdom and New York, 
NY, USA. 

Irish M, Irish G (2000) Agaves, Yuccas, and Related Plants. Timber Press, Portland, OR. 

Jensen ME, Burman RD, Allen RG (1990) Evapotranspiration and Irrigation Require- 
ments. American Society of Civil Engineers Manuals and Reports on Engineering 
Practices No. 70. ASCE, New York, NY. 

Kluge M, Schomburg M (1996) The Tonoplast as a Target of Temperature Effects in 
Crassulacean Acid Metabolism. In: Crassulacean Acid Metabolism, vol. 114 (eds 
Winter K, Smith JAC), pp. 72-77. Springer-Verlag, Berlin-Heidelberg. 

La Barre W (1938) Native American beers. American Anthropologist, 40, 224-234. 

Lewis S, Gross S, Visel A, Kelly M, Morrow W (2015) Fuzzy GIS-based multi-criteria 
evaluation for US Agave production as a bioenergy feedstock. GCB Bioenergy, 7, 
84-99. 

Li H, Foston MB, Kumar R et al. (2012) Chemical composition and characterization 
of cellulose for Agave as a fast-growing, drought-tolerant biofuels feedstock. RSC 
Advances, 2, 4951-4958. 

Mylsamy K, Rajendran I (2010) Investigation on physio-chemical and mechanical 
properties of raw and alkali-treated Agave americana fiber. Journal of Reinforced 
Plastics and Composites, 29, 2925-2935. 

NASS. (2015). National Agricultural Statistics Service, United States Department of 
Agriculture. Available at: http://www.nass.usda.gov/ (accessed 8 July 2015). 

Neales T (1973) The effect of night temperature on CO, assimilation, transpiration, 
and water use efficiency in Agave americana L. Australian Journal of Biological 
Science, 26, 705-714. 

Nobel PS (1991) Achievable productivities of certain CAM plants: basis for high val- 
ues compared with C3 and C4. New Phytologist, 119, 183-205. 

Nobel PS (1996) High Productivity of Certain Agronomic CAM Species. In: Crassu- 
lacean Acid Metabolism, vol. 114 (eds Winter K, Smith JAC), pp. 255-265. Springer- 
Verlag, Berlin-Heidelberg. 

Nobel PS (2003) Environmental Biology of Agave and Cacti. Cambridge University 
Press, Cambridge, UK and New York, NY, USA. 


Nobel PS, Smith SD (1983) High and low temperature tolerances and their relation- 
ships to distribution of agaves. Plant, Cell, and Environment, 6, 711-719. 

Nunez HM, Rodriguez LF, Khanna M (2011) Agave for tequila and biofuels: an eco- 
nomic assessment and potential opportunities. GCB Bioenergy, 3, 43-57. 

Osmond C (1978) Crassulacean Acid Metabolism: a curiosity in context. Annual 
Reviews in Plant Physiology, 29, 379-414. 

Owen N, Fahy K, Griffiths H (2015) Crassulacean acid metabolism (CAM) offers sus- 
tainable bioenergy production and resilience to climate change. GCB Bioenergy, in 
press doi: 10.1111/gcbb.12272. 

Pacheco-Sanchez C, Villa-Ayala P, Montes-Belmont R, Figueroa-Brito R, Jiménez- 
Pérez A (2012a) Effect of Ricinus communis extracts on weight and mortality of 
Scyphophorus acupunctatus (Coleoptera: Curculionidae). International Journal of 
Applied Science & Technology, 2, 83-94. 

Pacheco-Sánchez C, Villa-Ayala P, Montes-Belmont R, Figueroa-Brito R, Jiménez- 
Pérez A (2012b) Repellency of hydroethanolic extracts of Ricinus communis 
(Euphorbiaceae) to Scyphophorus acupuntatus (Coleoptera: Curculionidae) in the 
laboratory. Florida Entomologist, 95, 706-710. 

Post DF, Mack C, Camp PD, Sulliman AS. 1988. Mapping and characterization of 
the soils on the University of Arizona Maricopa Agricultural Center. Proc. Hydrol- 


ogy and Water Resources in Arizona and the Southwest, 49-60. University of Arizona, 
Tucson, Ariz. 

Saxton KE, Rawls WJ (2006) Soil water characteristic estimates by texture and 
organic matter for hydrologic solutions. Soil Science Society of America Journal, 70, 
1569-1578. 

Scanlon B, Reedy R, Stonestrom D, Prudic D, Dennehy K (2005) Impact of land use 
and land cover change on groundwater recharge and quality in the southwestern 
US. Global Change Biology, 11, 1577-1593. 

Scanlon B, Jolly I, Sophocleous M, Zhang L (2007) Global impacts of conversions 
from natural to agricultural ecosystems on water resources: quantity ver- 
sus quality. Water Resources Research, 43, W03437. doi:10.1029/2006WRO 
05486. 

U.S. Department of Agriculture (USDA) (2012) Pima Cotton 2012 Yield Per Har- 
vested Acre by County for Selected States. U.S. Department of Agriculture, 
National Agricultural Statistics Service. Available at: http://www.nass.usda.gov/ 
Charts_and_Maps/Crops_County/ (accessed 2 April 2014). 

van Velthuizen H, Huddleston B, Fischer G et al. (2007) Mapping Biophysical Factors 
That Influence Agricultural Production and Rural Vulnerability. Food and Agriculture 
Organization of the United Nations, Rome. ISBN 1684-8241. 

Winter K, Smith JAC (1996) Crassulacean Acid Metabolism: Current Status and Per- 
spectives. In: Crassulacean Acid Metabolism, vol. 114 (eds Winter K, Smith JAC), 
pp. 389-426. Springer-Verlag, Berlin-Heidelberg. 

Yan X, Tan DKY, Inderwildi OR, Smith JAC, King DA (2011) Life cycle energy and 
greenhouse gas analysis for agave-derived bioethanol. Energy & Environmental 
Science, 4, 3110-3121. 

Yang X, Cushman J, Borland A et al. (2015) A roadmap for research on crassulacean 
acid metabolism (CAM) to enhance sustainable food and bioenergy production in 
a hotter, drier world. New Phytologist, 207, 491-504. 


Supporting Information 


Additional Supporting Information may be found in the 
online version of this article: 
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